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Lattice-coupled antiferromagnetic spin model is analyzed for a number of frustrated lattices: 
triangular, Kagome, and pyrochlore. In triangular and Kagome lattices where ground state spins 
are locally ordered, the spin-lattice interaction does not lead to a static deformation of the lattice. 
In the pyrochlore structure, spin-lattice coupling supports a picture of the hexagon spin cluster 
proposed in the recent experiment [S. H. Lee et al, Nature, 418, 856 (2002)]. Through spin- lattice 
interaction a uniform contraction of the individual hexagons in the pyrochlore lattice can take 
place and reduce the exchange energy. Residual hexagon-hexagon interaction takes the form of a 
3-states Potts model where the preferred directions of the spin-loop directors for nearby hexagons 
are mutually orthogonal. 

PACS numbers: 75.10.Hk, 75.10.Jm 



The continued interest in insulating antiferromagnet 
appears to have a twofold objective. One is the search 
for a novel type of quantum ground state, other than 
the conventional Neel-ordered state, particularly in two 
dimensions^. A resonating- valence-bond ground state, if 
it should exist, is believed to lead naturally to a super- 
conducting phase when doped with holeA In a sepa- 
rate development, understanding the nature of antiferro- 
magnetic ground state defined on the frustrated lattice, 
even at the classical level, has been the focus of much 
theoretical^ as well as experimental activit}*^. 

In both these cases the fundamental Hamiltonian de- 
scribing the spin interaction is 



H 



(ii) 



Si ■ {Jij > 0) 
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defined for an appropriate set of bonds {ij). The an- 
tiferromagnetic interaction for an insulating magnet is 
mediated by the superexchange process in which over- 
lap of Wannier orbitals localized at different positions 
provides the necessary (virtual) hopping mechanism. As 
such, it is not surprising that the exchange energy Jij 
should depend on the separation of orbitals. Existing ex- 
perimental data suggests that Jij falls off as 6- 14th power 
of the separation^. 

Writing the equilibrium position of the ions by i and 
j, and the displacement vector of each ion by Ui and 
Uj^ the exchange integral has the expansion in the small 
displacements 

Jij = J{\i^Ui- j - Uj\) ^ Jo - JiCji • {uj - Ui) (2) 

where Cji = {j — — i\ is the unit vector, and Jo and 
Ji are positive constants. For an Einstein, or optical, 
phonon mode and ignoring the kinetic energy of the dis- 
placement, we arrive at the lattice-coupled spin model^ 
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^Y- 2 



(3) 



The purpose of this paper is to analyze the ground state 
of the model Hamiltonian, Eq. for a number of frus- 
trated lattices: triangular, Kagome, and pyrochlore^. 

By rescaling the displacements, Ui Uil\fK^ J\ is 
rescaled to J\l\fK = a, while the overall energy scale is 
fixed by Jo, which is set to one. The reduced Hamiltonian 
has the form 



H 
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Ot^Ui- fi^]-^U^^ 
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with fi 



Sj . We introduce the notation j G 



to indicate all bonds j that are exchange-coupled to site 
i. Minimizing the energy gives the condition relating the 
lattice positions with the spins: 



i/a={fi) = Y,eij{Si-S^). 



(5) 



Classical antiferromagnet on a triangular lattice has a 
ground state characterized by a 120° angle between a pair 
of adjacent spins. Long-range order is established at the 
mean-field level. When Eq. (|5|) is applied to the ground 
state spin arrangements of the triangular antiferromag- 
net, one finds Ui = {). Unrestricted numerical solution of 
Eq.(|5|) together with the classical mean- field equations 
for the spin average {Si) also yields = 0, consistent 
with a preliminary Monte Carlo simulation of the model 
Eq.(|3|!^. Thus, at the classical level, the ground state 
of the classical antiferromagnet on a triangular lattice 
is unaffected by the coupling to lattice. Employing the 
acoustic phonon model, {K/2) {ui — Uj)'^^ alters the 



.i{Ui 



^)/a = (fi), which 



mean- field equation to XljGiV 
also gives zero for u^^ . 

The effective Hamiltonian after "integrating out" the 
displacement Ui^ 



1 
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can be expanded around its classical minimum by use 
of the Holstein-Primakoff (HP) theory. Up to quadratic 
order in the HP bosons^^, 



Heff Eo 



k 



with 



7/c 
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E^ = --NzS\ 
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where 7fc = ^ E(i7 
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6 is the coordina- 



tion number for triangular lattice, and N is the number 
of sites. There are no terms in the effective Hamilto- 
nian proportional to up to quadratic order in the HP 
bosons, and the "vacuum" of the Hamiltonian remains 
unaltered after the introduction of spin-lattice coupling. 

Average of fi can be worked out within the quadratic 
theory, yielding 
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for the triangular lattice, and nB{k) is the boson occupa- 
tion number of energy cok- Therefore, small quantum or 
thermal fluctuation fails to produce a lattice distortion, 
or "spin-Peierls effects" as it is often known in the lit- 
erature. The higher-order terms in the HP bosons may 
have interesting consequences, including a spontaneous 
distortion of the lattice, and shall be considered in the 
future. 

Classical ground state of the antiferromagnet on the 
Kagome lattice is defined by the condition J^^^^ Si = 
where i = 1,2,3 are the three corners of a triangle r, for 
all r. This leads to the well-known ground-state degen- 
eracy of order 2^, where N is the number of triangles 
in the lattice. Nevertheless, ground-state spin-spin cor- 
relation 5*^ • Sj equals —1/2 for all nearest-neighbor pairs 
{ij). Although there is no long-range order, the spins are 
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FIG. 1: Ground-state spin configuration of the triangular and 
Kagome lattice. A,B,C refer to three orientations of the spin. 



locally ordered, with a coherence length of one lattice 
spacing. The spin ordering patterns for triangular and 
Kagome lattices are summarized in FigQ 

The r.h.s. of Eq.© gives zero for the Kagome lattice, 
as in the triangular lattice, because the local spin struc- 
tures are the same in both lattices. Si ■ Sj = —1/2. HP 
boson analysis of the Kagome lattice also reveals the ab- 
sence of static lattice distortion, whether with optical or 
acoustic phonons. Our findings may be summarized as 
follows: In the triangular lattice where spins are long- 
range ordered, and the Kagome lattice which has a local 
spin ordering without the long-range order, spin-lattice 
coupling fails to produce the static lattice distortion, or 
the spin-Peierls effect. 

Pyrochlore lattice is distinct from the previous two 
cases in the sense that spins do not order even locally. 
Ground state manifold of classical spins on the pyrochlore 
lattice is defined by the condition XliGr '^i — ^ where the 
basic building block r is a tetrahedron, for all the tetrahe- 
dra forming the lattice. The requirement is not sufficient 
to determine Si ■ Sj uniquely for nearest-neighbor sites, 
hence the claim that local spin ordering is missing in the 
pyrochlore lattice. 

Interaction energy — (a^/2) ff is minimized for ff 
maximum at each site. This selects the collinear spins 
(all spins parallel or antiparallel to each other) as the 
preferred ground state. The collinear spin patterns and 
the associated lattice distortion were analyzed in details 
in Ref. 0, following the pioneering work on the spin- 
Peierls effect in the pyrochlore structure in Ref. 0. 

Recently, neutron scattering data on the pyrochlore 
compound ZnCr204 revealed a very interesting picture 
of the spin dynamics^^. According to Ref. 0, each non- 
overlapping hexagon embedded inside the pyrochlore lat- 
tice has six spins form a collinear, antiferromagnetic clus- 
ter, which are energetically decoupled from those of other 
hexagonal clusters. The resulting block spins are chris- 
tened "spin- loop directors", or directors for short, in Ref. 

While this picture is intuitively appealing, no quan- 
titative justification for the formation, and the stability 
of, such a hexagonal spin cluster appears to exist to date. 

In the rest of the paper we show that spin-lattice cou- 
pling can aid in the formation of a hexagonal cluster sug- 
gested in Ref. 13. We argue that each hexagon can 
shrink - uniform contraction without breaking Cq sym- 
metry - to maximize the gain in exchange energy within 
the cluster. Spins within each hexagon is antiferromag- 
netically ordered, at the classical level. The contraction 
imposes a constraint, via Eq.Q, that the directors of 
nearby hexagons be mutually orthogonal, thus lifting the 
huge degeneracy of the original ground state manifold. 

We first analyze the case of a single hexagonal anti- 
ferromagnetic chain coupled to the lattice as in Eq.(|3|). 
The classical ground state of the effective Hamiltonian 
is given by the staggered spin, Sf = ±S'i. The stag- 
gered spins give Ui all pointing inward to the center of 
the hexagon, as given by 
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in the HP analysis. Here uj^ is the dispersion of a sin- 
gle hexagonal unit: uo^ = 2|sin/c|, k = 27r {integer)/ 6. 
The orientation of the antiferromagnetic spins defines the 
spin- loop director of Ref. 13. 

Taking individual hexagon as a structural unit, the 
pyrochlore lattice is built up of four different types of 
non-overlapping hexagons (defined A, B, C, and D), de- 
pending on the orientation of the face of a hexagon. 
The four directions correspond also to the normal of 
the four faces of a tetrahedron. Planes perpendicular 
to each orientation defines a Kagome lattice embedded 
in a pyrochlore structure. Taking one such Kagome 
plane as the basal plane, spanned by ei and 62, and 
the vector 63 connecting two nearby basal planes, the A- 
type hexagons are given the coordinates (3m,2n,p), or 
3mei + 2ne2 +^63, for integers m, n, and p (Fig|2|). The 
other three types are sandwiched between two nearby 
Kagome planes that are shown in Fig. [21 Coordi- 
nates of a B-type hexagon, for example, is defined by 
those of the A- type lying nearest to it. The remain- 
ing hexagons are then each located at {B; 3m, 2n + 
(C; 3m - 1, n, 2p + 1), and {D; 3m - 1, n, 2p). An A- 
type hexagon at (3m, 2n,p) is neighbored by four B-type 
hexagons located at (3m, 2n + l,p— 1), (3m, 2n — 1), 
(3m, 2n+l,p) and (3m, 2n— Two other hexagons, of 
type C and surround the A-type hexagon, too. Such 
"connectivity" of a given hexagon type to other hexagons 
can be worked out, for all hexagon types. 

A unit cell has a 2 x 2 x 3 structure with eight hexagons, 
two of each orientation^. Each hexagon is surrounded by 
six nearest-neighboring hexagons as shown in Fig. |31 Ac- 
cording to the previous discussion, a single hexagon un- 
dergoes a uniform contraction and has six collinear, an- 
tiferromagnetic spins. Requiring such hexagon contrac- 
tion throughout the whole pyrochlore lattice, the r.h.s. 
of Eq.Q imposes the condition that the directors of the 





FIG. 2: Three vectors that span the pyrochlore: ei, 62, and 
63. ei and 62 lie within a Kagome plane, while 63 connects 
the nearest-neighbor Kagome planes. Each plane is drawn by 
drawn by full and dashed lines. The angle between the three 
vectors are shown in a single tetrahedron. Thick black (first 
layer), and thick red (second layer) hexagons shrink due to 
spin-lattice coupling. 




FIG. 3: Local spin configurations of the distorted pyrochlore 
lattice. Short (long) bonds are denoted by color (dashed) lines 
after the lattice distortion. Each short bond belongs to one 
and only one hexagon cluster. 



nearest-neighbor hexagons be orthogonal (Fig|3|!"'^^. A 
self-consistency requirement on the hexagons is there- 
fore that the network of hexagons be "colored" in one 
of three colors, say R, G, and B, with no two neighbor- 
ing hexagons having the same color. In other words, the 
relative orientations of the directors satisfy an antiferro- 
magnetic 3-state Potts model 
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for the nearest-neighbour hexagons, (/J), and their re- 
spective directors, (dj^dj)^ where dj^dj takes on R, G, 
or B. We find, through explicit construction, that such 
coloring of the hexagon units can indeed be realized for 
a pyrochlore lattice. The strength of coupling between 
the directors Jg/ / depends on the spin-lattice interaction 
a, and also on the level of quantum fiuct nation within 
a hexagonal cluster. A severe quantum fiuct nation leads 
to the reduction of {Si ■ Sj) for z,j belonging to differ- 
ent hexagons, and an effectively weaker a in Eq.Q. In 
the extreme limit Jeff the hexagons are completely 
decoupled from each other, but more generally a resid- 
ual hexagon- hexagon interaction of order Jeff lifts the 
degeneracy and leads to a band of excitation spectra. 

So far the distortion mode is discussed in terms of the 
hexagon as a unit. Now, let us consider the distortion of 
an isolated tetrahedron. For a single tetrahedron, there 
are six vibrational modes: singlet Ai, a doublet and 
a triplet T^^^ . When the pyrochlore lattice assumes the 
hexagonal distortion we discuss, the resulting distortion 
for an isolated tetrhedron is not tetragonal {E) as dis- 
cussed in Ref. 11, but is a linear combination of a doublet 
E and a triplet T2 (FigI3|). 

In conclusion, lattice-coupled antiferromagnetic 
Heisenberg spin model on a variety of frustrated lattices 
is analyzed. For a triangular and Kagome lattice 
where the ground state spin structures are locally rigid. 
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FIG. 4: Distortion of an isolated tetrahedron consistent with 
a hexagon shrinkage. The spins configuration is coplanar. 

no spin-Peierls effect arises. Pyrochlore spins have a 
continuous manifold of ground states and a lack of local 
rigidity of spins, and spin-lattice coupling easily induces 
a lattice distortion leading to a reduction of the ground 
state degeneracy. A particular pattern of such lattice 
distortion, namely a hexagon contraction, is proposed 
and analyzed, following the suggestion of Refill. Each 
hexagonal cluster takes advantage of the lattice defor- 



mation to reduce its size and maximize the exchange 
energy within the cluster. The hexagonal ^protectorate^ 
of Ref 13 is thus obtained. Interaction between nearby 
hexagons is mediated by the spin-lattice condition, 
Eq.((n|), and leads to mutually orthogonal directors of 
nearby hexagons. In broader perspective, the spin-lattice 
interaction provides a channel for the self-organization 
of spins, which helps relieve the frustration inherent in 
the underlying microscopic Hamiltonian. 
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